.BLDIM is c(4,4), which results
in a 4 by 4 blocking dimension for distributed matrices.
Any time a function takes the bldim= argument, it will
default to this value unless the user specifies an alternative.
The default ICTXT is 0. This is the full 2-dimensional processor grid.